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Abstract. We study with semi-analytical methods the statistics of pronounced arcs caused 
by lensing of galaxies by foreground galaxy clusters. For the number density and redshift 
distribution of rich clusters we use Press-Schechter theory, normalized on the basis of empirical 
data. For the background sources we make use of observational results in the Hubble Deep 
Field. We present results for three different lens models, in particular for the universal profile 
suggested by Navarro, Frenk and White. Our primary concern is the dependence of the 
expected statistics on the cosmological parameters, Qm, ^Ia- The theoretical estimates are 
compared with the cluster arcs survey EMSS, and the resulting constraints in the f2-plane are 
presented. In spite of considerable theoretical an observational uncertainties a low-density 
universe is favored. Degeneracy curves for the optical depth and likelihood regions for the 
arc statistics in the Sl-plane depend only weakly on the cosmological constant. 
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1 Introduction 

There is now significant evidence that we are living in a critical universe in which 
the vacuum energy (or some effective equivalent) dominates and ordinary (baryonic) 
matter is only a tiny fraction of the total matter content: 

2 1 

The first two statements are strongly favored by the combination of recent measure- 
ments of the luminosity-redshift relation for type la supernovae at high redshifts and 
the observed temperature fluctuations of the cosmic microwave background radiation. 
(For a review, see 0.) JIa/ — 0.3 ± 0.1 is also consistent with data from rich clusters 
of galaxies, in particular from the analysis of the X-radiation emitted by the hot intra- 
cluster gas. The total matter contribution of rich clusters can also be determined by 
making use of (strong and weak) gravitational lensing. Such observations provide an 
estimate of the ratio ^Ib/^ai- a-nd Qm can then be obtained separately by using 
in addition the primordial abundances of light elements, synthesized during the first 
three minutes. (For a discussion of other methods, we refer again to and references 
therein.) 

Readers who do not closely follow these developments should, however, be warned 
that there remain serious worries. For instance, one cannot yet exclude systematic 
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effects in the supernovae data (e.g., intrinsic evolution, extinction) which could mas- 
querade a cosmological constant. 

In this situation any method which provides independent restrictions on the den- 
sity parameters ^Im, is of great interest. Lensing statistics on large scales has 
this potential. The use of gravitational lensing as a tool to determine cosmological 
parameters has been suggested in the pioneering works |^ . An early detailed study 
of lensing statistics was performed by Turner, Ostriker and Gott Much of the 
extensive later work is based on this paper. After some sporadic studies which investi- 
gated the effect of the cosmological constant on some specific aspects of lensing, it was 
pointed out in ||^, |^, |^ that statistical properties of gravitational lensing may exhibit 
strong A-dependencies. 

This is not the place to review the many more recent studies in the field of lensing 
statistics. We confine ourselves to the following remarks (more extensive references 
can be found in the quoted papers): 

(i) Several authors have recently re-analyzed the statistics of strong gravitational 
lensing of distant quasars by galaxies ^, 0. Observationally, there are only 
a few strongly lensed quasars among hundreds of objects. The resulting bounds on 

and Q\ are, however, not very tight because of systematic uncertainties in the 
galaxy luminosity functions, dark matter velocity dispersions, galaxy core radii and 
limitations of the observational material. 

(ii) On the basis of existing surveys, the statistics of strongly lensed radio sources 
has been studied in several recent papers Beside some advantages for 
constraining the cosmological model, there is the problem that the redshift distribution 
of the radio sources is largely unknown. (One can, however, make use of a strong 
correlation between the redshift and flux density distributions.) 

(iii) Clusters with redshifts in the interval 0.2 ^ ^ 0.4 are efficient lenses for 
background sources at Zs ~ 1. For several reasons one can expect that the probability 
for the formation of pronounced arcs is a sensitive function of ^Im and Q\. First, it 
is well-known that clusters form earlier in low density universes. Secondly, the proper 
volume per unit redshift is larger for low density universes and depends strongly on Q\ 
for large redshifts. An extensive numerical study of arc statistics has been performed 
by Bartelmann et al. |]l5|, with the result that the optical depth depends strongly on 
A. In our semi-analytical treatment, discussed below, we find, however, only a weak 
A-dependence. 

This paper is organized as follows. Since most relativists are not familiar with 
lensing statistics, we repeat in section |^ some basic concepts (cross sections, optical 
depth) and formulae needed in the theoretical analysis. For a comparison of theoret- 
ical expectations with observations, one has to include magnification effects and the 
observational selection criteria have to be taken into account. This will be explained 
in section ^ In section ^ we discuss further inputs we shall need in our semi-analytical 
treatment. This includes the Press-Schechter formula for the differential comoving 
cluster number density, dynamical and structural information required for the lens 
models, as well as redshift and luminosity distributions for the background sources. 
We shall compare our theoretical expectations with data for the Extended Medium- 
Sensitivity Survey (EMSS), which will briefiy be presented in section |^. Our main 
results are described in section g. These are compatible with a similar independent 
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investigation of Cooray jl^l , but do not agree with the claimed A-dependence of Bartel- 
mann et al. ||l5|. Possible sources of this conflict will be indicated. In section ^ we 
summarize and discuss future prospects for tighter constraints in the Q-plane. 



2 Generalities on lensing statistics 

As for any scattering situation, the concept of cross-sections is basic in statistical 
lensing. 

Consider a source at redshift Zs and distance Dg, as well as a lens a distance 
Dd (we always use angular diameter distances). We are interested in lensing events 
characterized by a set of lensing properties, Q, of this gravitational lens system. In this 
paper we shall consider especially lensing events in which a galaxy behind a cluster 
lens is seen as an arc with a length-to-width ratio R, larger than some given value 
Rq. For a given lens at distance Dd, the cross section Oq for the properties Q is by 
definition the area of the source sphere of radius , within which a source has to be 
located in order to be imaged with property Q. 

It is convenient to also introduce the corresponding cross section Oq in the lens 
plane by a simple rescaling: ctq = [Dd/ DsYgq. 

Consider next a distribution of lenses (deflectors) with number density nd{zd) at 
redshift Zd- The optical depth Tq for property Q is defined to be the probability that a 
given source at redshift z^ to undergo a lensing event with property Q. This is equal 
to the fraction of the entire source sphere of radius which is covered by the cross 
sections Uq (the cross-sections are in practice not overlapping): 



1 



nd{zd) 



dV{zd) 



dzd 



croizd) dzd- 



(1) 



Here, dV/dz denotes the proper volume per unit redshift. For a Fricdmann imiverse 
we have, in standard notation, 



with 



dy_ 

dz 

dt 

dz 



Ho {l + z)E{zy 



(2) 



(3) 



where 



E^{z) = Qm{1 + + {1- - f^A)(l + zf + ^A. 



(4) 



In connection with the geometrical factor (^) in the integral (Q) we recall that the 
angular diameter distance D{z) for a Friedmann universe is given by 



D(z) 



Ho{l + z) |n^|i/2 



dz' 

W) 



(5) 



where flu ■— I ^ fin ~ is a curvature parameter and S{x) = sin^, sinhx for 
if > 0, 0, < 0, respectively. The dependence of dV/dz on the cosmological parameters 
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is one of several sources for such a dependence of the optical depth. Others are implicit 
in the cross-sections and the lens distributions. Astrophysical uncertainties of these 
quantities limit attempts to exact cosmological parameters from lensing statistics. 

Inserting into (|l|), the optical depth can be expressed more simply in terms of 
the cross-section relative to the lens sphere: 

■^q(-Zs)= / nd{zd)<JQ{zd}\dt/dzd\dzd. (6) 
Jo 

Simple applications of this formula can be found in the early papers on lensing statis- 
tics. (See also jl^, especially chapters 11 and 12.) 

We now consider in more detail the statistics of arcs produced by cluster lenses. 
Let n{M, z) dM be the comoving number density of clusters at redshift z in the mass 
interval dM about M. The cross-section for the property that a background galaxy 
at redshift Zd behind a cluster of mass M at redshift z is seen as an arc with length- 
to- width ratio larger than i?o will be denoted by a{RQ,M,z,Zs). In terms of these 
quantities the optical depth (^ is given by 

t{>Ro,Zs)=( dz\dt/dz\{l + zf [ dMn{M,z)a{Ro,M,z,Zs). (7) 







Singular isothermal spheres (SIS) are the simplest lens models we shall use. For 
these the ratio R coincides with the magnification fi, and the cross-section is ll^ 



Here tT„ is the (constant) one-dimensional velocity dispersion, and Dds the angular 
diameter distance from the lens to the source. 



3 Amplification bias and selection 

For any magnitude limited sample of sources, the number of lensed sources is larger 
than it would be in an unbiased sample, because lensing brightens into the sample 
sources that would otherwise not be detected. This is, for example, a particularly pro- 
nounced effect in quasar lensing surveys, because the faint end of the quasar luminosity 
function rises steeply. 

Taking the amplification bias into account, the lensing probability (^ has to be 
replaced by 

Z"^" dr 

prob(> i?o, Zs,m-s) = / -7—B{ms, Zs, Zd) dzd, (9) 
Jo dzd 

where the first factor under the integral is the integrand in (|^) and the bias factor 
B, as a function of the apparent source magnitude ms and the redshifts Zs, Zd, can 
be expressed as follows: Let dNg/dms denote the unlensed number counts of sources, 
which agrees in practice with the observed number counts. Then 

/■°° dN 

B{ms,Zs, Zd) = {dNs/dms)-^ / -r^(m, -I- A, z,)g(A) dA, (10) 

. A__:__ dm. 
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where q{A) dA is the probabihty distribution for the magnification of the apparent 
magnitude of lensed images. Since A = 2.51og]^Q fi, the latter is simply related to the 
magnification probability distribution dfi. For a SIS model we have 

2 

= (^-1)3 ' Mmin = 2. (11) 

dNs/drris is simply related to the luminosity function of the source population, for 
which a Schechter form is often adopted. 

In optical arc searches only galaxies brighter than some limiting apparent magni- 
tude miim are counted. For a comparison of theory with observation we must include 
this observational selection in (||) and (|l^). Details of how this was implemented are 
described in |l^ . 



4 Additional theoretical and observational input 

For the numerical evaluation of (|^) and the corresponding giant arc statistics we need 
explicit knowledge of cross-sections, the differential comoving number density n(M, z) 
of clusters, and the distribution {dNs/dms){ms, Zg) of the source galaxies. These 
quantities are, unfortunately, not well-known. 

We have used three simple lens models, namely isothermal spheres with vanishing 
and positive core radii, and also the universal density profile proposed by Navarro, 
Frenk and White jl^. For isothermal models we need the velocity dispersion 
that enters in the cross-section with the fourth power (see (^)). ay depends on the 
cluster mass as well as on the cluster redshift z. For the ^-dependence we adopt the 
one of the spherical collapse model, while the M-dependence is derived from empirical 
correlations, using in particular the correlation between ay and the cluster temperature 
[ pO[ . (As expected from the virial theorem, the temperature is strongly correlated with 
the mass, T cx M^/^.) 



For n{M,z) we use the Press-Schechter formula |21 . This is based on a mixture 
of statistical reasoning, linear perturbation theory, and the spherical collapse model. 
Numerical simulations have demonstrated that the Press-Schechter expression provides 
a much better description than could be expected. One of the assumptions that goes 
into its derivation is that the smoothed density field is Gaussian. The function n{M, z) 
is essentially determined by the rms density fluctuations cr(M, z) of linear theory. (For 
textbook treatments, see |2^.) The z-dependence of a{M,z) is then equal to the 
one of the growing mode of linear perturbation theory, 

r 1 + z' 

Vgiz; Qm, - -nuEiz) J dz\ (12) 

and (7{M, z = 0) can be expressed in terms of the power spectrum of the density 
fluctuation in linear theory. For the latter we use the CDM transfer function given 
by Bardeen et al. ||2^, and assume that the initial fluctuations are adiabatic and 
nearly scale free. The normalization parameter erg is taken from psf . Uncertainties 
of this parameter, as well as its dependence on Qm, CIai are critical, since erg enters 
exponentially in the Press-Schechter formula. 
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For the source distribution of galaxies we make use of the present knowledge based 
on the Hubble Deep Field (HDF). The deep exposure of this survey is important for 
our analysis. Unfortunately, the redshift distribution is not precisely known, since a 
detailed spectroscopy of about 900 faint sources in the sample is impossible. There 
exist, however, two photometric redshift catalogs that were obtained by analyzing the 
observations in four broadband filters. We are always using the catalog in which 
is based on spectral template fitting. 

As a result of large K-corrections, the luminosity function is not well-known, and 
cannot be used for determining the amplification bias. Our method of treating this is 
discussed in 

The previous discussion should make it amply clear that we have to live with 
substantial uncertainties. Hopefully, these will be reduced in the near future. 



5 Extended Medium-Sensitivity Survey and arc statistics 

In the next section we shall compare our theoretical expectations with the results of a 
CCD imaging survey of gravitational lensing selected from the Einstein Observatory 
Extended Medium-Sensitivity Survey (EMSS). The latter is a large and sensitive X-ray 
catalog and has been studied by a number of groups, including followed up observations 
using the ROSAT or ASCA satellites. Presently, about 100 EMSS sources are classified 
as clusters. 

In a subsample was selected for an arc survey, using the University of Hawaii 
2.2 m telescope, and also the CFHT for specific observations. It should be emphasized 
that the EMSS is not strictly flux limited, since the detection of sources was optimized 
for high surface brightness. Therefore, more extended clusters with lower surface 
brightness, especially relatively close ones, may have been missed. Moreover, at high 
redshifts only the most luminous clusters can be seen. Since these are likely to be more 
massive, the lensing fraction is expected to become larger. For these reasons clusters 
in the redshift range 0.15 < z < 0.6, and surface brightness > A x lO^'' erg/s in 
the energy band 0.3 — 3.5 keV are most suitable for a comparison with theoretical 
expectations. 

The subsample selected in contains 21 such clusters. (This corresponds to 
about 1300 such clusters over the entire sky.) Among these there are 6 clusters with 
giant arcs satisfying the following criteria: (i) The length-to-width ratio R of the arc 
is > 10. (ii) The magnitude in the V-band satisfies mv < 22. (Actually, my was not 
always measured; in those cases < 22 was adopted, which should be about the 
same.) The 6 clusters in the limited area of the sample scale to about 360 over the 
entire sky, and the chance that a cluster in the specified range of z and contains 
an arc is roughly 0.3. 

So far the statistics are poor, but with the new satellites CHANDRA and XMM 
the situation should soon improve. 



R. Kaufmann et al., Giant Arc Statistics and Cosmological Parameters 



513 




Fig. 1 Differential optical depth for Zs = 1 and R > 10. Clusters are modeled as SIS. See 
text for details. 



6 Results 

In this section we present the main results of our study. A more detailed discussion 
and additional parameter studies are given in fl^ . 

Fig. 1^ shows the differential optical depths for the formation of giant arcs with 
i? > 10 for the following four cosmological models: SCDM (f^A/ = 1, f^A = 0), OCDM 
(flM = 0.3, f^A = 0), ACDM i^M = 0.3, = 0.7), rCDM {Qm = 1, f^A = 0, shape 
parameter T = 0.23). (The latter model mimics a scenario with an unstable massive 
r-neutrino.) The source is taken at redshift — 1. The foreground cluster lenses with 
X-ray luminosities in the EMSS band > 2 x 10*'^ erg/s, are modeled as singular 
isothermal spheres. 

Note - and this is one of our main points - that the curves for OCDM and ACDM 
do not differ much. The total optical depths are listed in the left column of Table 
|l|, w hile the right column contains those of the numerical study of Bartelmann et al. 
[[l5[ . We emphasize that the agreement is quite good for three of the four models, but 
for the OCDM model there is a difference of about an order of magnitude. Possible 
sources of this discrepancy will be discussed later. 

The dependence of the expected total number of giant arcs on the cosmological 
parameters 51^/ and J^a is shown in Fig. ^. Parameter variations and changes of the 
lens model have convinced us that the shape of the degeneracy curves (the lines of 
constant arc number) is quite stable, and we believe that their weak r^A-dependence 
will also be found in future more accurate studies. The absolute numbers refer to the 
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Fig. 2 Expected total number of giant arcs. > A x 10^^ erg s ^, niv < 22, ho = 0.7. 



EMSS subsample specified a the end of the previous section, and should be compared 
with the 360 arcs mentioned there. The best agreement is reached for fl^j « 0.4, 
practicahy independent of Oa. 

In Fig. 1^ we present the result of a maximum likelihood analysis, always keeping 
a SIS model. We have marginalized with respect to several empirical parameters, 
including as, which affect the prediction most strongly. At the present stage the 90% 
likelihood regions are still quite broad, but a low-density universe is clearly favored. 

We have not yet made a similar likelihood analysis for the universal profile of 



Table 1 Optical depth. Comparison between our results jl^ and the results in jl^. 



SCDM 3.1 X 10~* 4.4 x 10~* 

OCDM 2.9 X IQ-'^ 2.9 x lO"*^ 
ACDM 5.3 X 10"'^ 
1.6 X 10~^ 



15 



rCDM 



3.3 X 10"^ 

4.4 X 10"^ 
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0.8^ 




0.9 1 



Fig. 3 Likelihood contours. 



Navarro, Frenk and White ||l^. As already emphasized, the degeneracy curves are 
quite similar (the slopes become slightly negative), however, the absolute numbers turn 
out to be smaller |jl^. Other profiles have recently been proposed; their implications 
for the arc-statistics will be a subject of further studies. 

As already mentioned, Bartelmann et al. |l^ found in their N-body simulations 
for the optical depth a strong r^A-dependence (see Table |l|). It is, of course, true 
that a numerical approach provides in principle more realistic lensing properties of the 
clusters, since substructure and asymmetries are included. What we find, however, 
difhcult to understand is the claim in [|l5j that clusters are less concentrated for the 
ACDM model than for the OCDM model with the same Hm, and that this is the main 
reason for the lower number of arcs in the former case. A more recent simulation of 
another subgroup [ p8[ of the Virgo consortium indeed found no substantial differences 
of the cluster structures and mass profiles. 

We have investigated the various A-dependences which combine to the pattern of 
degeneracy lines shown in Fig. ^j. The geometrical factor in (|^) is, of course, common 
to all lensing statistics, and the A-dependence implicit in the SIS cross-section (^) is 
quite typical. For a conserved comoving lens number density, these would combine to 
degeneracy curves with JIa ~ ^m- (This kind of degeneracies are, for instance, found 
for the statistics of strongly lensed quasars.) This is changed by the A-dependence 
of the Press-Schechter formula for n{M, z) to the weak dependence in Fig. ^. (For 
further details, see JT8|.) 
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7 Conclusions and outlook 

Statistical strong lensing data for a specific class of sources have the potential to con- 
strain the cosmological parameters J7m and Oa- In combination with other information 
(from type la SNe, CMB anisotropics, etc.) this could tighten the allowed range of 
these important cosmological parameters, and also serve as a consistency check. 

So far, observational and theoretical uncertainties are limiting the use of lensing 
statistics for cosmology. Beside much better statistical samples we also need accurate 
knowledge about the population and mass distributions of the lenses, as well as on the 
redshift and luminosity distributions of the sources. 

In our semi-analytical study of giant arc statistics we have investigated the range 
of the current systematic uncertainties in constraining the cosmological parameters in 
the Um-^a plane. The available data for EMSS favor a low-density universe, but the 
90% upper bound for Qm is not really tight. More secure, we beheve, is our conclusion 
that the degeneracy curves in the fi-plane depend only weakly on A. In contrast to 
other claims in the literature , statistical lensing of giant arcs thus can not be used 
to constrain f^A- It niay, however, serve as a consistency check. (Note that the steep 
degeneracy curves are semi-transversal to those for type la SNe, as well as to the ones 
related to the position of the first acoustic peak in the CMB anisotropics.) 

One may be sceptical that strong lensing statistics will ever become an accurate 
tool. Weak lensing by large-scale structures looks more promising. The latter produce 
a tiny gravitational shear field which will be detectable in the near future. This shear 
pattern is correlated with the power spectrum of density fluctuations. The theoretical 
analysis is in this case more reliable than for strong lensing statistics, because linear 
perturbation theory is quite accurate on large scales. Weak lensing maps may soon be 
obtained with wide fleld imaging surveys. (For a review, and references, see p9| .) 

N. S. thanks the organizers of the Journees Relativistes 1999 for inviting him to beautiful 
Weimar and to present this work. He especially thanks G. Neugebauer for his gracious 
hospitality. N. S. thanks M. Bartelmann for a useful discussion on arc statistics. We are 
grateful to C. Frenk for providing us his FORTRAN subroutine for the calculation of the 
universal profile parameters. 
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